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Abstract 



In this paper we determine the velocity, the energy and estimate writhe and twist helicity 
contributions of vortex filaments in the shape of torus knots and unknots (toroidal and 
poloidal coils) in a perfect fluid. Calculations are performed by numerical integration of 
the Biot-Savart law. Vortex complexity is parametrized by the winding number w, given 
by the ratio of the number of meridian wraps to that of the longitudinal wraps. We find 
that for w < 1 vortex knots and toroidal coils move faster and carry more energy than a 
reference vortex ring of same size and circulation, whereas for w > 1 knots and poloidal 
coils have approximately same speed and energy of the reference vortex ring. Helicity 
is dominated by the writhe contribution. Finally, we confirm the stabilizing effect of 
the Biot-Savart law for all knots and unknots tested, that are found to be structurally 
stable over a distance of several diameters. Our results also apply to quantized vortices 
in superfluid 4 He. 

1 Introduction 

The study of vortex filament motion in an ideal fluid, and in particular of vortex rings in 
presence or absence of periodic displacements of the vortex axis from the circular shape (Kelvin 
waves), dates back to the late 1800s [39| 38J. Despite its long history, this study is still 
an active area of research [20[ [THl EE]. For example, the fact that upon large amplitude 
Kelvin waves vortex rings can slow down and even reverse their translational motion has been 
recognised only recently [221 [7]. Alongside the traditional interest for problems in classical fluid 
mechanics, additional interest is motivated by current work on superfluid helium [121 [9], [T71 [13] 
and atomic Bose-Einstein condensates [HI [23, [19] . 

Here we shall be concerned with slightly more complex vortex structures; namely vortex 
filaments in the shape of torus knots and unknots in an ideal fluid. Since these vortices 
are closely related to circular vortex filaments with small- amplitude distortions lying on a 
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mathematical torus, like vortex rings they propagate in the fluid by self-induction along the 
central axis of the torus, and they also rotate in their meridian plane (the poloidal plane of 
the torus) as their vortex core spins around the local center of mass, inducing an additional 
twisting motion of the vortex on itself, that cannot be neglected. 

Among all possible knot types, torus knots constitute a special family of knots amenable to 
particularly simple mathematical description. These knots can be described by closed curves 
wound on a mathematical torus a number of times in the longitudinal direction of the torus 
and a number of times in the meridian direction. If these two numbers are co-prime integers 
greater than one, then we have standard torus knots (see Figure H)), whereas if the curve 
winds the torus only once in one of the two directions, then the curve is simply unknotted 
(see Figure [2]), reducing to the standard circle when both numbers coincide with one. Even 
though torus unknots have trivial topology, their geometry may be rather complex, taking 
the shape of toroidal or poloidal coil, depending along which direction the curve is multiply 
wound. 

Vortex torus knots and unknots provide a good example of structures that are relatively 
complex in space, but still amenable to study relationships between dynamical properties, 
like velocity and energy, and geometric and topological features, a step further towards the 
study of more complex structures present in turbulent flows. The aim of this paper is thus 
to continue and extend previous work [35l 26J to more complex vortex structures. Here we 
shall investigate the propagation velocity and the kinetic energy of vortex torus knots and 
unknots in some generality, by comparing results to a standard vortex ring of same size. Since 
superfluid 4 He has zero viscosity, thus providing a realistic example of an Euler fluid, we shall 
choose circulation and vortex core radius as physically realistic quantities for quantized vortices 
in superfluid helium, and we shall carry out the research by direct numerical integration of 
the Biot-Savart law. Unlike previous vortex dynamics calculations of quantized vorticity 
[23 U El ED], however, we shall assume that no friction force [10J acts on the superfluid 
vortices; thus our results will apply to superfluid helium at temperatures below 1 K, where 
the dissipative effects of the normal fluid are truly negligible. 



2.1 Vortex motion under Biot-Savart and LI A law 

We consider vortex motion in an ideal, incompressible fluid, in an unbounded domain. The 
velocity field u = u(x, t), smooth function of the vector position x and time t, satisfies 



2 Mathematical background 



V • u = in IR 3 , 



u — 7- as x — > oo , 



(1) 



with vorticity uj defined by 



u = V x u 



V • u> = in IR 3 . 



(2) 
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In absence of viscosity, fluid evolution is governed by the Euler equations and vortex motion 
obeys Helmholtz's conservation laws [36]. Transport of vorticity is given by 

^ = V x (u x w) , (3) 

admitting formal solutions in terms of the Cauchy equations 

dx ■ 

u)i(x,t) = c^(a,0)^ . (4) 

From this expression we can see how both convection of vorticity from the initial position a 
to the final position x, and the simultaneous rotation and distortion of the vortex elements by 
the deformation tensor dxi/d&j are combined together. Since this tensor is associated with 
a continuous deformation of the vortex elements (by the diffeomorphism of the flow map), 
vorticity is thus mapped continuously from its initial configuration u>(a, 0) to the final state 
given by u>(x, £); hence the Cauchy equations establish a topological equivalence between initial 
and final configuration by preserving vorticity topology. In absence of dissipation, physical 
properties such as kinetic energy, helicity and momenta are therefore conserved along with 
topological quantities such as knot type, minimum crossing number and self-linking number 
[33]. 

The kinetic energy per unit density E is given by 

E = - I ||u|| 2 d 3 x = constant , (5) 
2 Jv 

where V — V is the fluid volume, and the kinetic helicity H by 



H — / u ■ u) d 3 x = constant . (6) 
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Here we assume to have only one vortex filament J 7 in isolation, where J 7 is centred on the 
curve C of equation X = X(s) (s being the arc-length of C). The filament axis C is given by 
a smooth (that is at least C 2 ), simple (i.e. without self- intersections), space curve C. The 
filament volume is given by V{J-) = 7ia 2 L, where L = L(C) is the total length of C and a is 
the radius of the vortex core, assumed to be uniformly circular all along C and much smaller 
than any length scale of interest in the flow (thin-filament approximation); this assumption 
is relevant (and particularly realistic) in the context of superfluid helium vortex dynamics, 
where typically a w 10 _8 cm. 

Vortex motion is governed by the Biot-Savart law (BS for short) given by 

. . r ri x (x-x(s)) , 

where T is the vortex circulation due to uj = Uot, where ujq is a constant and t = t(s) = dX/ds 
the unit tangent to C. Since the Biot-Savart integral is a global functional of vorticity and 
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Figure 1: Examples of torus knots with winding number w > 1 (left column) and w < 1 
(right column). Note that for given p and q, torus knot is topologically equivalent to the 
knot Tq, p , and we write T p , q ~ Tq tP - The tube centred on each knot is only for visualization 
purposes. 

geometry, analytical solutions in closed form other than the classical solutions associated with 
rectilinear, circular and helical geometry are very difficult to obtain. Considerable analytical 
progress, however, has been done by using the Localized Induction Approximation (LIA for 
short) law. This equation, first derived by Da Rios [H] and independently re-discovered by 
Arms & Hama |1] (see the review by Ricca [30] ) , is obtained by a Taylor's expansion of the 
Biot-Savart integrand from a point on C (see, for instance, the derivation by Batchelor [TT]). 
By neglecting the rotational component of the self-induced velocity (that in any case does not 
contribute to the displacement of the vortex filament in the fluid) and non-local terms, the 
LIA equation takes the simplified form 

Tc 

Ulia = — hi 5 b oc cb , (8) 

47T 

where c = c(s) is the local curvature of C, 5 is a measure of the aspect ratio of the vortex, 
given by the radius of curvature divided by the vortex core radius) and b = b(s) the unit 
binormal vector to C. 

2.2 Torus knots 

We consider a particular family of vortex configurations in the shape of torus knots in IR 3 . 
These are given when the curve C takes the shape of a torus knot 7~ p , q ({p, q} co-prime integers 
with p > 1 and q > 1), given by a closed curve wound on a mathematical torus p times 
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Figure 2: Examples of torus unknots with winding number w > 1 (poloidal coils; left column) 
and w < 1 (toroidal coils; right column). All these unknots are topologically equivalent to the 
standard circle. The tube centred on each unknot is shown only for visualization purposes. 

in the longitudinal (toroidal) direction and q times in the meridian (poloidal) direction (see 
Figure [1]) . When one of the integers is equal to one and the other is equal to m the curve is 
no longer knotted, thus forming an unknot homeomorphic to the standard circle Uq, but with 
a more complex geometry. Depending on which index takes the value m, the curve takes the 
shape of a toroidal coil (p = m) U m ^, or a poloidal coil (q = 1) lA\, m (see Figure [2]). When 
{p, q] are both rational, T PtQ is no longer a closed knot, the curve covering the toroidal surface 
completely. Here we shall consider only curves given by {p, q] integers. The ratio w = q/p 
denotes the winding number and Lk = pq the self-linking number, two topological invariants 
of T Pt q. Note that for given p and q the knot T p , q is topologically equivalent to T qtP , that 
is Tp,q ~ Tq, P i i.e. they are the same knot type, even though their geometry is completely 
different: toroidal coils resemble thick vortex rings with one full turn of twist of vorticity 
(internal swirl flow), whereas poloidal coils may recall circular jet flow generated by helicoidal 
distribution of vorticity on circular central axis. 

A useful measure of geometric complexity of C is given by the writhing number [TH] defined 

by 

Wr(C) = ±U i{8) X t(a * } • 1X(g) ~ X(S * )] ds ds* (9) 
{) -^%% ||X( S )-X( S *)||3 ' W 

where X(s) and X(s*) denote two points on the curve C for any pair {s, s*} G [0, L], the 
integration being performed twice on the same curve C. The writhing number provides a 
direct measure of coiling and distortion of the filament in space. This information can be 
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related to the total twist Tw and, by knowing the self-linking number Lk of each knot, to the 
helicity H of the vortex, that can be estimated to be H = TLk = T(Wr + Tw) [34], without 
resorting to direct calculation of the integral ([6]). 

Now, let us identify each knot T VA with the vortex filament (we shall hereafter refer to 
T p>q as the vortex torus knot), and consider dynamics and energy by using the BS ([7]) and 
the LIA law (jHJ) to determine the effects of different geometries and topologies on dynamical 
properties and kinetic energy. 

The existence of torus knot solutions to LIA were found by Kida [21] in terms of elliptic 
integrals. By re- writing LIA in cylindrical polar coordinates (r,a,z), and by using standard 
linear perturbation techniques, small-amplitude torus knot solutions (asymptotically equiva- 
lent to Kida's solutions) were derived by Ricca [28] . These latter give solution curves explicitly 
in terms of the arc-length s, given by 

r = r + e sm(w(f)) , 

a = — H cos{w(p) , , , 

r wr (10) 

t ( 1 \ 1/2 

z — he H 7T cos(u>0) , 

r V w J 

where r is the radius of the torus circular axis and e 1 is the inverse of the aspect ratio 
of the vortex. Since the LIA is related (30] to the one-dimensional Non-Linear Schrodinger 
Equation (NLSE), torus knot solutions ( 1T0|) correspond to helical travelling waves propagating 
along the filament axis, with wave speed k and phase <fi = (s — Kt)/ro- Vortex motion is given 
by a rigid body translation and rotation, with translation velocity u = z = u along the torus 
central axis and a uniform helical motion along the circular axis of the torus given by radial 
and angular velocity components r and a. In physical terms, these waves provide an efficient 
mechanism for the transport of kinetic energy and momenta throughout the fluid. 
By using eqs. (TTUT) . Ricca [25] proved the following linear stability result: 

Theorem 1 Let T VA he a small- amplitude vortex torus knot under LIA. Then T VA is steady 
and stable under linear perturbations if and only if q > p (w > 1). 

This result provides a criterium for LIA stability of vortex knots, and it can be easily 
extended to inspect stability of torus unknots (i.e. toroidal and poloidal coils). This stability 
result has been confirmed for the knot types tested by numerical experiments [35] • Interest- 
ingly, LIA unstable torus knot were found to be stable under the BS law, due to the global 
induction effects of the vortex. This unexpected result has motivated further work and some 
current research which is still in progress. 

3 Numerical method 

Dynamical quantities of vortex knots and unknots are evaluated by direct numerical integra- 
tion of the BS law (EJ) and direct numerical calculation of the other properties. The numerical 
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Figure 3: Examples of convergence test given by plots of the translational velocity u against 
the number of discretization points iV on the vortex axis. Similar results have been found for 
the other knots/unknots tested. Broken lines are for visualization purposes only. 



code is described in detail elsewhere [2] and it has been used also to study interaction and 
reconnection of vortex bundles [3]. The vortex axis is discretized into N segments and the 
Biot-Savart integral is de-singularized by application of a standard cut-off technique [371 HI E] • 
The time evolution is realized by using a 4 th order Runge-Kutta algorithm. Convergence has 
been tested in space and time as usual by modifying the number of discretization points and 
the size of the time step. 

In all calculations, the initial condition is given by ( FlUj) . We set r = 1 cm, e = 0.1, 
T = 1CT 3 cm 2 /s, (the value expected for superffuid 4 He), 5 = 2 x 10 8 /e 1 / 2 ^> 1 (typical of a 
superfluid helium vortex core radius a = 10~ 8 cm). To analyse unknots we replace in eqs. (fTUj) 
(1 — 1/w 2 ) with 1 1 — 1/ 1 . It is useful to compare these results with the dynamics of a vortex 
ring Uq of same size and vorticity. Thus we take a reference vortex ring with radius Tq — 1 cm. 
Typical time-step value in our calculations is 10~ 2 s. Convergence in time has been tested by 
ranging the time-step from 10~ 3 s to 5 x 10~ 2 s. Examples of convergence in space are shown 
in Figure [3] for torus knot 72,5 an d unknot U^\. The calculations are performed using constant 
mesh density N/ L chosen following convergence tests similar to those shown in Figure |3J We 
set N/L = 50 for poloidal unknots U\ im (m = 2, 3,..., 7), N/L = 20 for toroidal unknots 
Um,i ( m — 2, 3, . . . , 7), N/L = 30 for knots 72, g (q = 3, 5, 7, 9) and N/L = 7 for knots T p ^ 
(p = 3,5,7,9). Typical errors in computing the velocity and the energy are approximately 
10~ 5 cm/s and 10 -7 cm 5 /s 2 , respectively. For the reference vortex ring we set the winding 
number w — 1 and iV = 313; then L = 6.26 cm, thus N/L = 50, translational velocity 
u = 1.38 x 10~ 3 cm/s, kinetic energy E = 8.77 x 10~ 6 cm 5 /s 2 and obviously H = cm 4 /s 2 . 
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Figure 4: Vortex knots: numerical values of the calculated quantities. Entries in the last 
column report whether (Y) or not (N) a reconnection event has taken place on or before time 
t. 
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Figure 5: Vortex unknots: numerical values of the calculated quantities. Entries in the last 
column report whether (Y) or not (N) a reconnection event has taken place on or before time 
t. 
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4 Results: velocity, energy and helicity 



Numerical values of calculated quantities are reported in the tables of Figure H] and El Let 
us consider first purely geometric information, that will be useful to understand the following 
results on velocity, energy and helicity. 

4.1 Total length and writhing number 

Diagrams of total length and writhing of knots and unknots are shown in Figure EJ From 
diagrams of total length (left column), we see that these plots reflect the elementary fact 
that longitudinal wraps contribute to total length more than meridian wraps. The marked 
difference in the slope of the two plots in each diagram is just due to the dominant contribution 
to L by the longitudinal wraps compared to the modest one of the meridian wraps. This 
behavior will be reflected in the relative kinetic energy of the system. 

Similarly for the amount of coiling and distortion of the filament axis in relation to the 
dynamics and helicity of the vortex. As we can see from the diagrams on the right of Figure [6j 
meridian wraps contribute modestly (if not at all appreciably to our degree of accuracy) to the 
total writhing of the filament, the dominant contribution coming from the longitudinal wraps 
present. Moreover, since the self-linking of the vortex filament is a topologically conserved 
quantity [31] given by Lk = Wr + Tw, where Tw denotes total twist of the filament, then 
information on writhing provides direct information on total twist, by taking Tw = Lk — Wr. 
Furthermore, being the number of meridian wraps directly proportional to total twist, we 
have an immediate estimate of the relative contributions to helicity (see the discussion in the 
sub-section below). 

4.2 Translation velocity 

The translational velocity u = z along the central axis of knots and unknots is calculated by 
using the Biot-Savart law. Absolute values are reported in the tables of Figure H] and [51 The 
diagrams of Figure [7] show the normalized velocity u/uq of knots and unknots plotted against 
the winding number. The velocity is greatly influenced by the relative number of longitudinal 
wraps, which contribute greatly to the total curvature of the vortex. In general the velocity 
decreases with increasing winding number. Fastest vortex systems are thus torus knots with 
highest number of longitudinal wraps. In the case of unknots we can see that meridian wraps 
actually slow down poloidal coils Wi m , making them traveling slower than the corresponding 
vortex ring. At very high winding number, torus knots and poloidal coils seem to reverse 
their velocity, thus traveling backward in space. This curious phenomenon has been observed 
independently [22j 13 > an d it can be justified on theoretical grounds by information based on 
structural complexity analysis [HI [32] . 

An estimate of the relationship between normalized velocity u/uq and winding number w 
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Figure 6: Total length L and writhing number Wr plotted against the winding number w 
of the torus knots and unknots calculated (p = 3,5,7,9; q = 3,5,7,9; m = 2, 3,..., 7). 
The isolated circle denotes the reference vortex ring value. Broken lines are for visualization 
purposes only. 

of torus knots is obtained by a linear regression, given by: 

f 4.41 -9.13w + 7.09w 2 , (w<l), 

U/Un = < . (11) 

1 \ 1.25 - 0.04w O > 1) , v ' 

with standard deviation of 0.031. 

In the limit w — > 0, the knot covers the toroidal surface completely with an infinite number 
of longitudinal wraps, and vorticity becomes a sheet of toroidal vorticity, with the induced 
velocity, purely poloidal in the interior and exterior region of the torus, that jumps across the 
sheet in opposite directions. The regression (TIT]) suggests a theoretical limit value u/u = 4.41 
that seems to be independent of the aspect ratio of the torus. In the other limit, w — > oo, the 
knot covers the torus with an infinite number of meridian wraps, and in this case vorticity 
induces a toroidal jet in the interior region. 
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Figure 7: Translational velocity u normalized with respect to the velocity Uq of a reference 
vortex ring plotted against the winding number w of the torus knots and unknots calculated 
(p = 3, 5, 7, 9; q = 3, 5, 7, 9; m = 2, 3, . . . , 7). The isolated circle denotes the reference vortex 
ring value. Broken lines are for visualization purposes only. 

4.3 Kinetic energy 

The diagrams of the normalized kinetic energy per unit density E/E (Figure El left column) 
of the knots and unknots calculated by using the Biot-Savart law, and the normalized energy 
calculated by the LIA law Eua/E (Figure M, right column), plotted against the winding 
number w, show trends dictated by the behavior of total length L = L(w) (cf. the diagrams 
of Figure El). Here the calculation of the volume integral (J5j) is replaced by the more economical 
line integral [8] 

E=- iu-Xxtds. (12) 
2 Jc 
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Figure 8: Left column: kinetic energy E calculated by using the Biot-Savart law and normal- 
ized with respect to the energy E of the reference vortex ring, plotted against the winding 
number w of the knots and unknots calculated (p = 3, 5, 7, 9; q = 3, 5, 7, 9; m = 2, 3, . . . , 7). 
Right column: kinetic energy -Exia calculated by using the LIA law and normalized with re- 
spect to the energy E plotted against w. The isolated circle denotes the reference vortex ring 
value. Broken lines are for visualization purposes only. 



An estimate of the relationship between the normalized energy E/Eq of torus knots and their 
winding number w is given by a non-linear regression based on least squares: 

E/E = 2.35 + 67.03e~ 6 - 44w , (w < 1) , (13) 

with root mean square deviation of 0.218. 

Different trends and values are obtained by calculating the normalized energy by using the 
LIA law; in this case, by using (jSJ), we have 

£LIA = i/^||u LIA || 2 d»X=(^) 2 ^^d S , (14) 

that is one of the conserved quantities associated with the LIA law [28J. Direct comparison 
between the diagrams of the two energies reveals two distinct trends: for w < 1 the LIA law 
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Figure 9: Translational velocity u of some knots and related unknots reported here for com- 
parison. In general knots travel faster than their corresponding unknots. The isolated circle 
denotes the reference vortex ring value. Broken lines are for visualization purposes only. 



under-estimates the actual energy of the vortex (knotted or unknotted), whereas for w > 1 the 
LIA provides much higher energy values. The importance of non-local effects is evident: the 
LIA energy of 7g,2 and i/^i, for example, is about 40% less than the corresponding BS energy, 
whereas 72,9 and Uij under LIA have 3 and 14 times more energy than their corresponding 
BS counterparts. These differences are essentially due to the contributions from the induction 
effects of nearby strands, captured by the BS law, but completely neglected under LIA. 

Another quantity of interest is the kinetic energy per unit length E/L. From the tables of 
Figures H] and [5] we see that for w < 1 E/L can vary by up to 100%. This result has interesting 
implications for the interpretation of experiments on superffuid turbulence, particularly at 
low temperatures, a regime that is dominated by Kelvin waves [24] . What is measured in the 
experiments is the total vortex length per unit volume A, and the turbulence energy per unit 
volume is deduced by multiplying the observed vortex length per volume times the energy 
per unit length; by integrating the velocity field of a straight vortex filament, the energy per 
unit length is estimated as E/L ~ r 2 /(47r 2 ) In (b/a), where b m A -1 / 2 is the typical distance 
between vortices. Our results show that for highly bent vortex filaments (as in the case of 
superfluid turbulence, where vorticity is even fractal [23]), E/L is certainly not constant. 



4.4 Helicity 

It is interesting to investigate the effects of topology on the velocity and kinetic helicity, by 
exploring the interplay of geometric and topological aspects on the dynamics. For this let us 
consider Figure [HI where the velocities of a torus knot and that of the corresponding unknot 
with same number of longitudinal or meridian wraps are shown for comparison. As we see 
in general vortex knots travel faster than their corresponding unknots, and the higher the 
number of longitudinal wraps the faster is the motion, whereas the higher the number of 
meridian wraps the slower is the propagation speed. 

As far as helicity is concerned, we have H = TLk [34] and for torus knots we can set 
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Figure 10: Plots of space travelled (z) and time elapsed (£) before a reconnection event takes 
place for torus knots T Pt 2 (p = 3, 5, 7, 9) plotted against the winding number w. Calculations 
are based on the Biot-Savart law. Broken lines are for visualization purposes only. 



Lk = pq, that is invariant for each knot type (Tw being the total twist) and it increases with 
knot complexity. Since Lk = Wr + Tw from the Wr- values of the table of Figure H] we can 
compute Tw for each knot. In general twist helicity is larger for knots with higher number 
of meridian wraps. For the unknots we can set Lk = Wr + Tw = (for the reference vortex 
ring we have zero writhe and zero twist); thus from the values of the table of Figure we 
have Tw = —Wr and therefore we can estimate writhe and twist helicity contributions. 

4.5 Structural stability considerations 

Aspects of structural stability of vortex knots and unknots based on permanence of knot 
signature and occurrence of a reconnection event are explored. In previous work we noticed 
the stabilizing effect of the BS law on LIA- unstable knots [32] • Here we extend this comparison 
to the knots and unknots considered so far. In general a vortex structure is said to be stable 
if it evolves under signature-preserving flows that conserve topology, geometric signature and 
vortex coherency. Let us consider results reported in the tables of Figure H] and [5] and shown 
in Figure [LTU 

Since our results concern both Euler's and superfluid dynamics, it is important to re- 
mark that superfluid vortices can reconnect with each other in absence of dissipation [251 02] , 
whereas in the Euler context vortex topology is conserved; this important difference has been 
recently reviewed by Barenghi [5]. Here we adopt the following criterium: when, during the 
evolution, a vortex reconnection takes place, then we stop the calculation, and deem the struc- 
ture to be unstable. The last column of the table of Figure @] and refers to the occurrence of 
a reconnection at instant t. In the absence of reconnection, we report the distance z travelled 
by the vortex during the computational time t. If z is much larger than the typical vortex 
size, the knot/unknot tested is said to be structurally stable. 
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From the last two columns of Figure H] we see that for w < 1 the space travelled tends to 
decrease with increasing knots complexity; for example, til^^) 3> tilg^)- The effect is also 
shown in Figure [101 Note that these knots would be LIA- unstable. Thus, the stabilizing effect 
due to the BS law is confirmed: these knots can indeed travel a distance which is larger than 
their size before unfolding and reconnecting. 

5 Conclusions 

In this paper we examined the effect of several geometric and topological aspects on the 
dynamics and energetics of vortex torus knots and unknots (i.e. toroidal and poloidal coils). 
This study is carried out by numerically integrating the Biot-Savart law, and by comparing 
results for several knots and unknots at different winding numbers w. Generic behaviors are 
found for the class of knots/unknots tested, and main results are presented by normalizing 
velocity and energy by the corresponding value of a standard vortex ring (Uq) of same size 
and circulation. 

In general, for w < 1 (where the number of longitudinal wraps is larger than that of 
meridian wraps) the more complex the vortex structure is, the faster it moves, and both torus 
knots and toroidal coils move faster than Uq. For w > 1 (where the relative number of meridian 
wraps dominates) all vortex structures move essentially as fast as Uq, almost independently 
from their total twist. Therefore, for all the structures tested total twist provides only a 
second-order effect on the dynamics. 

We have also found that for w < 1 vortex structures carry more kinetic energy than Uq, 
whereas for w > 1 knots and poloidal coils have almost the same energy as Uq. The LIA law (an 
approximation often used to replace the Biot-Savart law, more computationally demanding) 
under-estimates the energy of knots with w < 1 and over-estimate the energy for w > 1. 

Kinetic helicity, that is naturally decomposed in writhe and twist contributions, is evidently 
determined by the relative number of longitudinal and meridian wraps present, the latter 
contributing to twist helicity rather modestly (relatively). 

Finally, we can extend previous results and confirm that for w < 1 the Biot-Savart law 
has a stabilizing effect on knots that are LIA-unstable: all vortex structures tested have been 
found to be structurally stable regardless of the value of w, being able to travel in the fluid 
for several diameters before eventually unfolding and reconnecting. 
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